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ABSTRACT 

We present the surface brightness profile of M31's stehar halo out to a projected radius of 175 kpc. 
The surface brightness estimates are based on confirmed samples of M31 red giant branch stars derived 
from Keck/DEIMOS spectroscopic observations. A set of empirical spectroscopic and photometric 
M31 membership diagnostics is used to identify and reject foreground and background contaminants. 
This enables us to trace the stellar halo of M31 to larger projected distances and fainter surface 
brightnesses than previous photometric studies. The surface brightness profile of M31's halo follows 
a power-law with index —2.2 ± 0.2 and extends to a projected distance of at least ^ 175 kpc 2/3 
of M31's virial radius), with no evidence of a downward break at large radii. The best-fit elliptical 
isophotes have b/a = 0.94 with the major axis of the halo aligned along the minor axis of M31's 
disk, consistent with a prolate halo, although the data are also consistent with M31's halo having 
spherical symmetry. The fact that tidal debris features are kinematically cold is used to identify 
substructure in the spectroscopic fields out to projected radii of 90 kpc, and investigate the effect 
of this substructure on the surface brightness profile. The scatter in the surface brightness profile is 
reduced when kinematically identified tidal debris features in M31 are statistically subtracted; the 
remaining profile indicates that a comparatively diffuse stellar component to M31's stellar halo exists 
to large distances. Beyond 90 kpc, kinematically cold tidal debris features can not be identified due 
to small number statistics; nevertheless, the significant field-to-field variation in surface brightness 
beyond 90 kpc suggests that the outermost region of M31's halo is also comprised to a significant 
degree of stars stripped from accreted objects. 

Subject headings: galaxies: halo — galaxies: individual (M31) — galaxies: structure 



1. INTRODUCTION 

Stellar halos of galaxies are low density environ- 
ments in which the detritus of hierarchical struc- 
ture formation can remain visible for Gigayears in 
the form of tidal debris structures. Detailed sim- 
ulations of stellar halo format io n in a cosmologi- 
cal c o ntext (iBullock et all 120011: iBuUock fc JohnstonI 
20051; IPont et al.' '2006^, '20081 'Johnston et al.l 12008!; 
Zolotov ct al. 2009, 2010; Cooper et al. 2010; Font et al.] 
201 1( ) are providing a growing framework for interpret- 
ing observations of tidal debris in stellar halos. Current 
simulations suggest that galaxy mergers play a primary 
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role in the formation of stellar halos: the outer regions of 
stellar halos are built via the accretion of smaller galax- 
ies and tidal stripping of their stars, while the inner re- 
gions may be built through a com bination of rnerging 
events and in s itu star formation (" Zolotov et al.l 120101 : 
iFont et al.ll2011 : McCarthy ct al. 2012^. Thus, the global 
properties of stellar halos are a product of the merging 
history of the host galaxy. 

Recent observational campaigns have extended 
studi es of st el lar h al os beyond the Local Group 
(Tan aka et al.l 1201 1[ iMartinez-Delgado et aD 120101 : 
Rad burn-Smith et al. 120111 ). However, M31 remains 
one of the best laboratories for studying stellar halos. 
M31's proximity enables us to study its resolved stellar 
population not just with photometry, but also through 
spectroscopy of individual stars. 

Early work on t he integrate d light of M31 

(|de VaucouleursI 119581 : [mhcrbos fc KennicuttI Il987| l 
was followed by resolved stellar population studies 
over the last several deca des (e.g., iMould fc Kristiaiil 
Il986t iFerguson et al.l [200^ . These studies found that 
within i? ~ 30 kpc, the M31 stellar halo appeared 
to be a continuation of M31's inner bulge, and was 
characterized by the following properties: a de Vau- 
couleurs r^/'^ law radial surf ace brightness profile 
(jPritchct fc van den Berghl 119941 ). a high stellar density 
(IReit zel et al. I 119981 ). and a high mean metallicity 
([Fe/Hl ^ -0.5 Mould & Kristian 1986; Rich ct 
iigMlDmrcll et al.ll200ll: iRertzel fc Guhathakurtal 12001 
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iBellazzini et al.ll2003t fPurrell et al.ll2004[ ). Furthermore, 
deep imaging (obtained with the Advanced Camera for 
Surveys on the Hubble Space Telescope) along M31's 
minor axis at radii of 10-35 kpc (well beyond the 
extent of MSl's disk) reve aled a sig nificant population 
of ^ 6-8 Gyr stars (.Brown et alJ [200l [2006, 200Z, 
|2008() . The properties of the inner regions of MSl's 
spheroid stand in stark contrast to the properties of 
the Milky Way's (MW) stellar halo, which consists 
almost entirely of old, metal-poor stars and appears to 
follow a power-law density profile {p{r) oc r'^) . Studies 
of the MW's halo have typically measured 7 ^ — 2.7 
to -3.5 (e g- iMorrisonet al l [20001: l Yannv et all [2OOOI: 
iSiegel et alJ l2002t lJuric'eta l 2008; Ses ar et all l20Tl[ ). 
equivalent to a surface density distribution with a slope 
of -1.7 to -2.5. 

In the past five years, large photometric and spectro- 
scopic surveys have discovered and begun to characterize 
MSl's outer stellar halo, which has a power-law surface 
brightness profile, low stellar density, an d metal-poor 
stars. Using the spectroscop ic data set of [Gilbert et al.l 
(|2006( ). iGuhathak urta et al.l (|2005f ) showed that beyond 
projected radial distances of R 



proj 



20-30 kpc, the sur- 
face brightness profile of M31's stellar halo is consistent 
with a power-law of index —2.3 out to the limits of the 
surveyed region, i?proj 165 kpc. A concurrent photo- 
metric study based on data from the Isaac N ewton Tele- 
scope Wide-Field Camera (jlrwin et al.l 120051 ) also found 
evidence for a break in the radial surface brightness pro- 
file of M31 at i?proj '~ 20 kpc, and determined that be- 
yond this radius the profile fiattens and is consistent 
with a power-law index of —2.3; their study reached to 
a proj ected radial distan ce of 55 kpc from the center of 
M31. Ilbata et al.l (|2007[ ) presented a large photometric 
survey of M31's southern quadrant undertaken with the 
Canada-France-Hawaii Telescope/MegaCam and found 
that the radial surface brightness profile of M31 out to 



-1.9. 



130 kpc is consistent with a power-law of index 



M31's outer stellar halo consists of stars that are on av- 
erage more metal-poor than the stars that com prise the 
inner, bulge-li ke spheroid. Kalirai et al.l ()2006af ) used the 
data set from iGilbert et al.l (|2006D to analyze the metal- 
licity of red giant branch (RGB) stars from i?proj ~ 11 
to 165 kpc in MSl's southern quadrant. They found 
that while the inner regions of MSl's spheroid are rel- 
atively metal-rich, the stellar population becomes in- 
creasingly m etal- poor beyond Rprr,\ ^ SO kpc (see also 
iKoch et al.ll2008l: iTanaka et al.ll2010l) . A contemporane- 
ous study found evidence of metal-poor MSI RGB stars 
from i?pioj ~ 10 to 70 kpc in fiel ds located near MSl's 
major axis ([Chapman et al.ll2006( ). 

Most previous studies of MSl's stellar halo consisted 
primarily of fields in the southern quadrant of MSI or 
along MSl's southern minor axis. Two recent photomet- 
ric surveys have extended observations of MSl's stellar 
halo into the other three quadrants, providing a more 
global view of the halo. The Pan- Androm eda Archaeo- 
logical Survey (PAndAS; iMcConnachie et al. 2009) used 
the Canada-France-Hawaii Telescope (CFHT) and the 
MegaCam instrument to survey the eastern, northern, 
and western quadrants of MSl's halo out to a maxi- 
mum projected radial distance of ^ 150 kpc. PAn- 
dAs found stellar sources consistent in color and ap- 



parent magnitude with RGB stars at the distance of 
MSI th roughout the survey volume ([McConnachie et all 
I2009D . ITanaka et all (|2010D surveyed along MSl's mi- 
nor axis with the Suprime-Cam instrument on the Sub- 
aru Telescope, achieving a deeper photometric limit than 
PAndAS, although over a much more limited area. They 
observed strips of contiguous fields on the southeast 
and northwest minor axis out to projected distances of 
-Rproj ~ 60 and ~ 100 kpc, respectively, and found the 
radial surface brightness profile along MSl's northwest 
minor axis is described by a power-law with an index of 
—2.2, consistent with results from the southeast minor 
axis. 

Interpretation of the global properties of MSl's stel- 
lar halo is complicated by the presence of large. 



extended, 



tidal 



a nd n u merous 
Ibata et al.l 12001': 'Fer guson et al 



2006bt Ilbata et al. 2007 



d ebris features 
[20021: 



; atures (e.g., 

S Kalirai et aLl 
120071 I2009r 



McConnachie et al.l [20091: rTanaka et al.ll2010[) . It is dif 



ficult to assess the effect of substructure on measured 
values (mean metallicity, surface brightness profile) with 
purely photometric studies. The surface brightness pro- 
file of the smooth underlying stellar population must be 
determined by removing entire areas with photometri- 
cally identified substructure. In much of MSl's stellar 
halo, this requires removing a significant fra ction of the 
available spatia l area from the analysis (^e.g.. Ilbata et al.l 
2007; Tanaka et al.l [20Tol ) . Spectroscopy of the stellar 
population yields the kinematical information needed 
to separate intact, dynamically cold tidal debris fea- 
tures (typically with ay ^ 30 km s~^) from an un- 
derlying spatially diffuse, dynamically hot population 
{ay > 100 km s^^). 

This paper is the first in a series that describes the 
global properties of MSl's stellar halo using spectro- 
scopic observations of member RGB stars in all quad- 
rants of MSl's stellar halo. In this contribution, we mea- 
sure the surface brightness profile of MSl's stellar halo 
using counts of spectroscopically confirmed MSI RGB 
stars and quantify the effect of substructure on the sur- 
face brightness profile. The profile presented here encom- 
passes data from more than 100 spectroscopic 4'xl6' slit 
masks in 38 fields; 27 fields are being presented for the 
first time. The fields cover a range of position angles and 
projected distances from MSI, and were all obtained with 
the DEIMOS spectrograph on the Keck telescope. The 
data included in our analysis are discussed in Section [2j 
The method used to convert counts of spectroscopically 
identified MSI RGB stars to a surface brightness estimate 
is described in Section [3| The global surface brightness 
profile of MSl's stellar halo is presented and analyzed in 
Section [4l The results are discussed within the broader 
context of hierarchical stellar halo formation in Section[5l 
Finally, we summarize our conclusions in Section [6] A 
distance modulus of 24.47 is assumed for all conversions 
of angular to physi cal units (correspon ding to a distance 
to MSI of 783 kpc: lStanek fc Garnavich 1998) . 

2. DATA 

We present data in 38 fields spanning a large range 
in position angle and projected distance from the center 
of MSI (Figure [H and Tabled]). Thirteen fields lie out- 
side the southern quadrant (the most heavily observed 
quadrant of MSl's halo), and 10 fields lie at projected 
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distances of 100 kpc or more from M31's center. The 
data were obtained as part of the SPLASH (Spectro- 
scopic and Photometric Landscape of Andromeda's Stel- 
lar Halo) survey of M31's stellar halo. In this section, 
we summarize the photometric and spectroscopic obser- 
vations and data reduction (Section [23J and the identi- 
fication of the M31 RGB star sample (Sections 12.21 and 
12. 3|) . Further details of the observations and reduction 
techniques used in the SPLASH survey can be found in 
the listed references. 

2.1. Photometric and Spectroscopic Observations 

The photometric data are primarily from imaging ob- 
servations obtained with the Mosaic Camera on the Kitt 
Peak National Observatory (KPNO) 4-m Mayall tele- 
scoped Photometric catalogs were derived from obser- 
vations in the Washington system M and T2 filters an d 
the intermediate-width DD051 filter (jOstheimed 120031 ). 
The observed magnit udes were transformed to V and / 
using the equations of iMaiewski et all (|2000[ ). 

Additional photometric data in the outer (i?proj > 
30 kpc) fields came from two sources. The photometry 
in field 'dlO' was derived from V an d / images obtaine d 
with the William Herschel Telescope (|Zucker et al.ll2007f ). 
Photometry for fields 'streamE' and 'streamF' was de- 
rived from V and / images obtained with the Suprime- 
Cam instrument on the Subaru Telescope (|Tanaka et al.l 
[MS). 

Photometric catalogs for the innermost fields (i?proj < 
30 kpc) were derived from observations obtained with 
the MegaCam instrument on the 3.6-ni CFHTF^ Images 
were obtained with the g' and i' filters. The observed 
stellar magnitudes were transformed to Johnson-Cousins 
V and / magnitudes using observations of Landolt pho- 
tometr ic standard stars as described in iKalirai et al.l 
(l2006aD . 

The above photometric catalogs were used to design 
spectroscopic slit masks for use in the DEIMOS spec- 
trograph on the Keck II telescope. Table [T] lists the 
filters used for each of the spectroscopic fields. Ob- 
jects were assigned a priority for inclusion on the slit 
mask based on image morphology, their Iq magnitude, 
and (when available) the probability of their being RGB 
stars at the distance of M31 based on their position in an 
(M- T2, M-DD051) color -color diagram (jPalma et al.l 
120031: IMaiewski et al.ll2005D. The mask design pro cess 
is described in detail bv iGuhathakurta eFaH (l2006h . In 
the high-density inner regions of M31's spheroid not all 
high priority M31 RGB candidates within the DEIMOS 
mask area can be included due to slit conflicts. In con- 
trast, the low density outer fields contain a high fraction 
of low priority filler targets due to the paucity of M31 
RGB candidates. The effect of this radial dependence of 
slit mask targets on the surface brightness estimates is 
discussed in Section r3. 2. 21 

Kitt Peak National Observatory of the National Optical As- 
tronomy Observatory is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

MegaPrime/MegaCam is a joint project of CFHT and 
CEA/DAPNIA, at the Canada- France-Hawaii Telescope which is 
operated by the National Research Council of Canada, the Insti- 
tut National des Science de I'Univers of the Centre National de la 
Recherche Scientifique of France, and the University of Hawaii. 



Spectra were obtained over nine observing seasons 
(Fan 2002-2010) with the DEIMOS spectrograph on the 
Keck II 10-m telescope. The 1200 line mm~^ grating, 
which has a dispersion of 0.33 A pixel" ^, was used for all 
observations. The slit width of 1" yields a resolution of 
1.6 A FWHM. The typical wavelength range of the spec- 
tral observations is 6450-9150A. This wavelength range 
includes the Can triplet absorption feature at ^ 8500 A 
and the Nai absorption feature at 8190A. Masks were 
observed for approximately 1 hr each, with modifications 
made for particularly good or bad observing conditions. 

The spectra were reduced using modified versions 
of the spec2d and specld softw are developed at the 
University of Calif ornia, Berkeley ([Newman eFalllml 
iCooper et al.ll2012l ). The spec2d routine is used to per- 
form the flat-fielding, night-sky emission line removal 
and extraction of one-dimensional spectra from the two- 
dimensional spectral data. The specld routine cross- 
correlates the resulting one-dimensional spectra with 
template spectra to determine the redshift of the ob- 
jects; the template library includes stellar spectra ob- 
tained with DEIMOS and galaxy spectra f rom the Sloan 
Digital Sky Survey (see Isimon fc Gehal 120071 for de- 
tails). Each spectrum was visually inspected, and only 
stellar spectra with secure velocity measurements wer e 
included in the analyzed data set ([Gilbert et al.ll2006[ ). 
Two corrections were applied to the observed velocities: 
(1) a heliocentric correction, and (2) a correction for im- 
perfect centering of the star within the slit. The latter 
correction was calculated by measuring the observed po- 
sition of the atmospheric A-band absorption feature rel- 
ative to night sk y emission lines ([Simon fc Gehal I2OO7I : 
ISohn et al.ll2007l) . 

We are presenting data from 108 spectroscopic masks. 
These masks tareted ~ 11000 objects, and yielded suc- 
cessful velocity measurements from over 5800 stellar 
spectra (52.2%). The remainder of the targets were 
galaxy spectra (20.2%), spectra for which a velocity 
measurement was not possible due to insufficient signal- 
to-noise ratio (S/N) or a lack of strong spectral lines 
(24.2%), and a small fraction of catastrophic failures 
(3.4%). 

2.2. Classification of M31 RGB Stars and Foreground 
MW Stars 

MW dwarf stars along the line of sight constitute a 
significant fraction of the observed stellar spectra even 
though the sample of spectroscopic targets consists pri- 
marily of objects that have been photometrically prese- 
lected to be M31 RGB candidates. Moreover, the line- 
of-sight velocity distributions of M31 halo RGB and MW 
dwarf stars overlap, and the velocity distribution of fore- 
ground MW stars has a tail that extends well into the 
velocity range typical of M31 halo stars. To estimate 
the surface brightness of a given M31 field, the relative 
fractions of stars in each of the two populations must be 
securely measured. 

Individual stars are identified as M31 red giants or 
MW d warfs using the method developed bv lGilbert et al.l 
([20061) . to which readers are referred for full details of 
the technique. A combination of photometric and spec- 
troscopic diagnostics are used to determine the proba- 
bility an individual star is an M31 red giant or MW 
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Fig. 1. — Locations and orientations of the spectroscopic masks (red rectangles) included in the analysis (Section [21l • The footprints of 
the KPNO/Mosaic (black), CFHT/MegaCam (dark gray), and Subaru/Suprime-Cam (light gray) images used to design the spectroscopic 
slit masks are also shown. The center of M31 is marked with an open circle, and the orientations of the major and minor axes of M31's 
disk are denoted respectively by the long and short solid lines. The dwarf elliptical satellites of M31 are denoted by circles, and the dwarf 
spheroidal (dSph) satellites of M31 are denoted by triangles. The dotted circles denote radii of 2, 4, 6, 8, and 11 deg from M31's center. 



dwarf: (1) line-of-sight velocity (wios), (2) photometric 
probability of being a red giant based on location in 
the (M — T2, M— DD051) color-color diagram (when 
available; Table [J), (3) the equivalent width of the Nai 
absorption line (surface-gravity and temperature sensi- 
tive) versus {V — T)o color, (4) position in the (/, V — I) 
CMD and (5) spectroscopic (based on the EW of the 
Ca II triplet) versus photometric (comparison to theoret- 
ical RGB isochrones) [Fe/H] estimates. Each diagnostic 
provides separation between M31 RGB stars and MW 
dwarf stars based on different physical parameters. A 
star's location in each diagnostic is compared to empir- 
ical probability distribution functions (PDFs) to deter- 
mine the likelihood {Li — logiQ(FRGB/^'dwarf)) the star 
is an M31 red giant or MW dwarf (Figure [5]); each PDF 
was determined us ing training sets of M31 red giant and 
MW dwarf stars (jGilbert et al.1 [20061) . The likelihoods 



for each diagnostic are combined to give the overall like- 
lihood, (Li), the star is an M31 red giant or MW dwarf. 
Stars that are > 3 times more probable to be an M31 red 
giant than MW dwarf star {(Li) > 0.5) are designated as 
secure M31 red giants (and vice versa for secure MW 
dwarf stars: (Li) < —0.5); stars below this probability 
ratio threshold are designated as marginal M31 red gi- 
ants (0 <{L,)< 0.5) or MW dwarfs (-0.5 <{Li) < 0). 
The analysis is restricted to stars designated as secure 
M31 red giants or MW dwarfs because th e secure cate- 
gorie s suffer from the least contamination (jGilbert et alj 
I2OOI . 

Figure[2]presents the position in the four most powerful 
diagnostics of secure M31 RGB and MW dwarf stars in 
fields interior (left panels) and exterior (right panels) to 
^proj ~ 80 kpc. Very few stars are identified as secure 
M31 RGB stars in the outermost fields, however Figure[2] 
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Fig. 2. — Locatio ns of stars in four of the five empirical diagnostics fSection l2.2IIGilbert et al.ll2006l ): (a) line-of-sight velocity, (b) DD051 
parameter (Section [2TTJ, (c) Nal EW, and (d) spectroscopic versus photometric [Fe/H] estimates. Left 4 panels: stars identified as secure 
M31 red giants (classes +2 and +3) and MW dwarfs (classes —2 and —3) in fields from Rproj ~ 8 to 80 kpc. The majority of the stars 
come from the interior fields (-Rproj < 30 kpc). Right 4 panels: stars identified as secure M31 red giants and MW dwarfs in the outer 
fields (-Rproj ^ 80 kpc). The RGB and dwarf probability distribution functions shown in panels (a) and (b) are normalized to equal area. 
The PDF contours in panels (c) and (d) show the 90%, 50%, and 10% contours. Although significantly fewer in number and primarily 
metal-poor (as can be seen by comparing panel (d) in the left and right figures), stars identified as M31 RGB stars in the outer fields are 
consistent with the distribution of M31 RGB stars in the inner fields. 



shows that those that are identified as M3f red giants 
are more consistent with the M3f RGB, rather than the 
MW dwarf, PDFs. The RGB stars in the outer fields 
(i?proj ^ 80 kpc) are bluer (more metal-poor) than the 
average RGB star in the inner fields, as expected based 
on the observed meta l hcity gradient of MSl' s stellar halo 
(iKalirai et alJl2006al : iChapman erall l2006t iKoch et alJ 
I2008D . However, their locations in the four diagnostics 
are consistent with the distribution of M31 RGB stars in 
the inner fields. 

The noticeable asymmetry in the velocity distribution 
of M31 RGB stars within 80 kpc (Figure [i;a), left panel) 
is due to two effects. The first is the presence of stars 
associated with the giant southern stream in many of 
the south quadrant M31 halo fields. The velocity of the 
giant southern stream varies with i?proj, but remains 
more negative than M31's systemic velocity, v^^^ = 
-300 km s~^ (llbata et al.l 120051 iGuhathakurta et alJ 
I20M IKalirai et al.ll2006bl : IGilbert et al. I12QO90. The sec- 
ond effect is incompleteness in the secure sample of 
M31 RGB stars at the positive end of the M31 veloc- 
ity distribution. The velocity diagnostic typically pro- 
vides greater separation of the populations, and thus has 
intrinsically higher weight, than the other diagnostics. 
Thus the velocity region where the M31 and MW veloc- 
ity distributions overlap tends to have a higher percent- 
age of stars identified as marginal M31 RGB s tars (for a 
full discussion, see the appendix of IGilbert et al. 2007). 

2.2.1. Comparison of M31 RGB and MW Dwarf Star 

Samples 



The left panels of Figure [3] show the distribution of 
stars in likelihood space. As discussed above, stars with 
marginal M31/MW classification (shaded region of left 
panels) tend to lie in the region of velocity space where 
the RGB and MW dwarf star PDFs overlap. Although 
the sample is divided into M31 and MW stars at {Li)=0, 
the true likelihood distributions of M31 {{Li)> 0) and 
MW stars {{Li)< 0) are expected to have tails that ex- 
tend into negative and positive likelihood values, respec- 
tively. Restricting the sample to secure M31 {{Li)> 0.5) 
and secure MW {{Li)< —0.5) stars minimizes the con- 
tamination in the M31 and MW samples. 

The right panels of Figure [3] display the location of 
stars in the plane defined by heliocentric line-of-sight 
velocity and photometric metalhcity ([Fe/H]). [Fe/H] 
values are determined by comparison of each star's /q 
magnitude and [V — I)q color with a grid of theoreti- 
cal RGB isochrones at the distance of M31 (Figure EJ. 
For stars that are either slightly bluer than the bluest 
isochrone, or slightly above the TRGB, the metallicities 
are measured by extrapolating the grid, enabling a con- 
tinuous mapping of the color-magnitude distribution of 
stars to [Fe/H]phot space. Clearly, these isochrone based 
[Fe/H] phot estimates are only physically meaningful for 
the bona fide M31 RGB stars. Regardless, in the [Fe/H]- 
vios plane, three distinct groups of stars are apparent in 
the sample. As discussed below, these groups represent 
RGB stars in the halo of M31, dwarf stars in the MW 
disk, and main-sequence turn off stars in the MW halo. 

The M31 RGB star locus is centered at v\os — 
—300 km s^^ (the systemic velocity of M31) and spans 
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Fig. 3. — Left: likelihood distribution fSection I2.2| l of all observed stars as a function of projected distance from M31's center. Blue 
histograms show the distribution of stars classified as M31 RGB stars, while black histograms denote foreground MW dwarf stars. The 
shaded area shows the region of likelihood space excluded from the analysis: stars that have likelihood values in this range are considered 
marginal identifications. Some stars are classified as MW dwarf stars even though {Li) > 0. These are stars whose Vq, Iq magnitudes place 
them bluer than the most metal-poor isochrone by an amount greater than the photometric error (leading to an extrapolated [Fe/H]pijot< 
—3.5, see right panels). Right: distribution of stars in the [Fe/H]pjiot~'yios plane. Marginal stars are those with —0.5 <(ii) < 0.5. Although 
the relative fraction of M31 RGB stars and MW dwarf stars changes with radius, the locus of M31 RGB stars in the [Fe/HJphot^fios plane 
remains populated even at large projected distances (-Rproj > 100 kpc) from M31's center. M31's systemic velocity is vi^c^ = —300 km s^^. 



a reasonable range of [Fe/H]phot for RGB stars. Even 
in the most distant fields (i?proj > 100 kpc), the sample 
includes stars that are clearly located in the M31 locus. 

The two MW dwarf star loci denote the two primary 
populations of MW contaminants expected in the sur- 
vey. The MW dwarf locus at velocities — 150 < vios < 
50 km s~^ and relatively red colors (which translate into 
high [Fe/H]phot estimates) is consistent with the proper- 
ties expected for stars belonging to the disk of the MW. 

The second MW dwarf locus has a wide spread in ve- 
locities (—400 < "yios S 0) ^iid significantly bluer colors 
than either the M31 RGB stars or MW disk stars; in- 
deed, the {V — I)o colors of these stars place them well 
to the blue of the RGB fiducials, and thus these stars 



are depicted with upper limits. This locus is consistent 
with the properties expected for main sequence turn-off 
stars at varying line-of-sight distances in the Galactic 
halo. The stars with {Li)> in Figure |3] that are clas- 
sified as MW stars are part of this population. These 
stars' velocities (wios ^ —200 km s~^), combined with 
their blue colors [{V — /)o^ 1], typically lead to high 
{{L.,)> 0) likelihood values because the CMD, Nai EW, 
and [Fe/H] diagnostics have no discriminating power at 
these colors (Figure [2]). Therefore, any star that is blue- 
ward of the most metal-poor isochrone by more than the 
typical photometric error is classifie d as an MW d warf 
star (Figured see Section 4.1.2 of .Gilbert et al.|[2006i. for 
a more detailed discussion). It is possible that a small 
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Fig. 4. — Io,{V — I)o color— magnitude dia gram for the spec- 
troscopic sample, with theoretical isochrones of IVandenBerg et al.l 
((2006) ranging from [Fe/H]=-2.3 (left most) to 0.49 (right most) 
overlaid (assuming an age of 10 Gyr and [a/Fe]=0.0). The thick 
isochrones denote metallicities of [Fe/H]=— 1.0 and 0.0, while the 
dashed isochrone is for an age of 5 Gyr and [Fe/H]=— 2.3. Stars 
classified as red giants at the distance of M31 trace the RGB 
isochrones, while stars classified as MW dwarf stars along the line of 
sight have a broader color distribution. Typical photometric errors 
as a function of / magnitude are shown on the left: stars that are 
bluer than the most metal-poor isochrone by more than the typica l 
photometric error are classified as MW dwarf stars (Section l2.2.1|l . 

fraction of these could be M31 stars. Individual stars 
might have spurious photometric measurements, or they 
could be asymptotic giant branch stars or very metal- 
poor M31 halo red giants. However, the vast majority 
of these stars do appear to lie well-removed from the lo- 
cus of M31 stars and to be part of a separate population 
of stars, with properties consistent with main sequence 
turn-off stars in the Galactic halo (Figure [3]). 

2.3. Identification of M31 Field Stars in Masks 
Targeting dSph Galaxies 

The SPLASH survey includes an effort to measure 
the kinematical and chemical properties of M31's dwarf 
satellite galaxies u s ing spectroscopy of memb er stars 
(|Kalirai et al.l [20091 [2010l : iTollerud et alllMl) . Thus, 
a number of the M31 spectroscopic fields presented here 
were designed to target dSphs in the halo of M31 (Ta- 
ble [1]). The number of field stars belonging to M31's 
stellar halo in the dS ph masks is e s timat ed following 
the technique used by [Gilbert et al.l ()2009l ) (their Sec- 
tion 2.5.1, Figures 3 and 4). The distribution of stellar 
objects in velocity, metallicity, and sky position is used to 
isolate RGB stars that are likely part of M31's stellar halo 
rather than members of the dwarf satellite. To classify 
the stars, we leverage the following three properties of the 
dSph galaxies: (1) they are generally compact enough to 
cover only a portion of the spectroscopic slit mask, (2) 
they have small velocity dispersions (< 10 km s"'^), and 
(3) their member stars generally span a limited range of 
[Fe/H]. All RGB stars bey ond the es timated King lim- 
iting radius of the dSph ([McConna chie & Irwin 2006; 



stars. RGB stars that are inside the limiting radius of 
the dSph but well removed from the tight locus in v\os - 
[Fe/H] space occupied by the dSph members are also 
designated as M31 field stars. Table [1] lists literature 
references for each dSph, where the distribution of stars 
and their membership is discussed in detail. 

Although this procedure will miss M31 field stars that 
happen to fall both within the limiting radius of the dSph 
galaxy and in the region of wios~[Fe/H] space occupied 
by the dSph members, the number of such M31 field 
stars is expected to be small. This is due in part to 
the fact that in each field, the dSph members span a 
very narrow range in velocity space compared to M31 
halo stars. It is also due to the fact that the surface 
brightness of each dSph becomes increasingly dominant 
over the surface brightness of M31's spheroid toward the 
center of the dSph, resulting in an increased likelihood 
that spectroscopic slits will be placed on a dSph member 
rather than an M31 field star. 

3. SURFACE BRIGHTNESS ESTIMATES FROM 
M31 RGB STAR COUNTS 

The surface brightness of M31's stellar halo is esti- 
mated from counts of spectroscopically confirmed M31 
RGB stars (Section 12. 2p . The ratio of securely identi- 
fied RGB stars (A^msi) to securely identified MW dwarf 
stars (A^Mw) in a field is used to estimate M31's stel- 
lar surface density. However, the density of foreground 
MW dwarf stars is not constant over the widely spaced 
spectroscopic fields. To account for this, the ratio of ob- 
served M31 RGB and MW dwarf stars is multiplied by 
the total expected surface density of MW stars in each 
field ba sed on the Besanc on Galactic population model 
fA^RnM^ lRobin et al]|2003h . The Besangon model counts 
are computecF^ for the same magnitude and color range 
as the spectroscopic targets, and standard model param- 
eters were used (models were drawn from a 1 deg^ line- 
of-sight centered on each mask, including all ages and 
spectral types). 

The scaled RGB counts are also corrected for two ob- 
servational sources of bias: pre-selection of likely M31 
RGB candidates using M, T2, and DD051 photometry 
(cddo; Equation (|3]), Section l3.2.ip . and unequal sam- 
pling of the RGB luminosity function (clf! Equation ([5|), 
Section I3.2.2p . Finally, the scaled and corrected RGB 
counts are converted to an 7-band surface brightness. 
The normalization is determined by fitting for an over- 
all normalization factor, /norm, that best matches our 
surface brightness estimates to the min or axis J-band 
surfac e brightness profile publishe d by Courteau et al.l 
(|2011l ) (dis cussed further in Section 2T|). The normaliza- 
tion of the ICourteau et al. | profile is based on an /-band 
M31 image obtained bv IChoi et al.l ([200^. The surface 
brightness estimates are thus calculated as 
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(1) 



Martin et al.' '2006'; 'M ajewski et al.l l2007i: iZucker et al.l 
2007; .Collins et al.. 2Qm are designated as M31 field 



Our method of estimating the surface brightness in 
each field is based on counting the relative numbers of 
M31 RGB and MW dwarf stars. Selecting only securely 
identified M31 RGB and MW dwarf stars (Section [22]) 

http://model.obs-besancon.fr/ 
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Fig. 5. — Examples of multi- Gaussian fits (solid curves) to the 
velocity distribution of M31 stars (bold histogram) in fields with 
kinematically cold tidal debris. These fits are used to determine 
the fract ion o f stars belonging to M31's kinematically hot spheroid 
(Section I3.1|l . Thin histograms denote the velocity distribution 
of all stars with recovered velocities. The fields shown span a 
range of projected distances from M31's center and demonstrate 
the range of velocity dispersions measured for the kinematically 
cold components detected in individual fields (dotted curves) , from 
~ 5 km s~^ (e.g., panel (d)) to ~ 50 km s~^ (panel (a)). In some 
fields, more than one kinematically cold component is identified 
(panel (c)). In all fits, the mean velocity and dispersion of the 
Gaussian representing M31's kinematically hot spheroid (dashed 
curves) is held fixed with parameters (i))''P'^= —300 km s~^ and 
ctI^^= 129 km s-l JGilbert et al.ll2007l) . 

is one way of doing this counting. We have also tested 
including the marginal M31 RGB and MW dwarf can- 
didates, and have investigated fitting the M31 RGB and 
MW dwarf star likelihood distributions (Figure [3]) to find 
the fraction of stars in each population. These alter- 
nate counting methods produce similar results. In the 
end however, our choice was influenced by the fact that 
restricting the sample to stars securely identified as be- 
longing to one population or the other provides a clean 
sample from which to fit for kinematically cold compo- 
nents in each field. 

3.1. Multiple Stellar Populations 

Many of the M31 halo fields targeted during the course 
of our Keck/DEIMOS survey show evidence of multi- 
ple kinematical components within the M31 halo popu- 
lation: i.e., one or more distinct, kinematically cold com- 
ponents associated with tidal debris in ad dition to M31's 
kinematically broad spheroid ([Gu hathakurta et al.ll20d6l : 
IKalirai et al.ll2nn6bl : IGilbert et al...2nn7. ,20091) . The ef- 
fect of these fields on the surface brightness profile of 
M31's stellar halo will be discussed in detail below. 

In fields with kinematically cold components (Table [1]), 
the number of M31 stars in the dynamically hot compo- 
nent is calculated by multiplying the total number of 
M31 RGB stars in each field by the fraction of stars esti- 
mated to belong to the kinematically hot component. 
This fraction is determined by calculating maximum- 
likelihood multi-Gaussian fits to the stellar velocity dis- 
tribution. The velocity distributions and associated 



maximum-likelihood fits of previously published fields 
with substructure are described in detail in the litera- 
ture (Table [J). 

Figure [5] shows examples of maximum- likelihood Gaus- 
sian fits in fields spanning a range of projected distances 
from M3I's center. The analysis implicitly assumes the 
presence of a dynamically hot component. In all fields 
with identified kinematical substructure, an underlying 
distribution of stars with a wide range of velocities is 
observed. However, this does not preclude the possibil- 
ity that the broad distribution of stars (particularly in 
fields at larger radii) is a superposition of fainter, in- 
distinct tidal debris features, rather than a dynamically 
well-mixed population of stars. 

Only two of the new M31 halo fields have kinemat- 
ical evidence of tidal debris. The velocity distribu- 
tion of stars in these fields is presented in Figures [SJ^c) 
and (d). The substructure in field d9 is at (u)''"'^= 
-329.7t6°6^ km s^^, with al''^ ^12. 6]}^% km s^^; 

29.ltlti% of stars in this field are estimated to be- 
long to the kinematically cold component. This field 
is located at 37 kpc near M31's major axis. The ex- 
pecte d disk velocit y at this location is ~ —100 km s~^ 
tibat a et al.l |2005[) . making it unlikely that this fea- 
ture is related to M31's disk. The peak of MW stars 
with V ~ —130 km s~^ may be related to the TriAnd 
overdensit y in the MW halo; t his feature is discussed 
further in iTollerud et al.l ()2012D . The substructure in 
field R06A220 is at (?;)"^''= -373.5 ± 3.0 km s-\ with 
cr«"b^g j+2.7 g-i. 42%±i5% of stars in this field are 
estimated to belong to the kinematically cold component. 

Finally, we note that of the 39 fields presented here, 
only 4 specifically targeted tidal debris features: f207 
and H13s targeted the giant southern stream, while the 
streamE and streamF fields w ere chosen to t arget photo- 
metric overdensities noted in iTanaka et al.l ((2010. ) . The 
spectroscopic data show clear kinematical substructure 
in f207 and II13s, but there is no detectable kinematical 
substructure in the streamE or streamF fields. The re- 
maining 35 fields were observed without prior knowledge 
of the presence of substructure in these fields, and thus 
represent a random sampling of the properties of M31's 
stellar halo. 

3.2. Primary Sources of Bias and Systematic Error 

Since the surface brightness estimates are based on 
counts of spectroscopically confirmed samples of M31 
RGB and MW dwarf stars, nonuniform spectroscopic tar- 
get selection and variable observing conditions can sig- 
nificantly bias our measurements. Radially dependent 
systematic errors that affect the slope of the measured 
surface brightness profile are the most important to con- 
sider, since a purely empirical normalization is used to 
convert counts of confirmed M31 RGB stars and MW 
dwarf stars to surface brightness estimates. The sur- 
face brightness estimates are corrected for two sources of 
bias: selection of likely M31 RGB stars for spectroscopy 
based on Washington AI, T2, and DD051 photometry, 
and the luminosity function (LF) of spectroscopic targets 
with successful velocity measurements. We describe each 
of these in detail below (Section 13.2.11 and sec:lumfun). 
Other sources of error that are not accounted for are dis- 
cussed in Section [3731 
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3.2.1. DD051 Selection Ef, 



iciency 



Selection of likely M31 RGB stars using photometry in 
the Washington M and T2 filters and the DD051 fiher 
(Section l2.1l) boosts the measured RGB/dwarf star ratios 
relative to fields without DD051 photometry, and results 
in an offset in the surface brightness between fiel ds with 
and without DD051 photometry ( Guhathakurt a et al.l 
120051 ). Among the masks that benefit from DD051-based 
target selection, masks at large i?proj contain a higher 
fraction of "filler" targets (that fail the DD051 criterion 
for RGB stars) than masks at small i?proj • Since many of 
these filler targets are MW dwarfs, this effect reduces the 
measured M31 RGB/MW dwarf ratio. Uncorrected, this 
would cause the derived radial surface brightness profile 
to be steeper (a more negative power-law index) than the 
true profile because a larger fraction of slits will target 
likely M31 RGB stars in DD051-based masks at small 

^proj • 

To correct for this effect, we estimate the factor by 
which the measured RGB /dwarf star ratios are increased 
by DD051-based target selection. For each field a 
DD051 selection function (DSF) is computed, defined 
as the ratio of the number of stars selected for inclusion 
on a spectroscopic mask (iVtarget) to the number of stars 
available (A'totai), as a function of the DD051 parameter 
(/ddo5i; Figure El^b)). The /ddo51 value measures the 
probability a star is an M31 red giant based on the Wash- 
ington M, T2, and DDO Sl photometry (jPalma et al.l 
[200I iMaiewski et al.|[2005l ). The DD051 correction fac- 
tor (cddo) is then defined as the ratio of two integrals: 
(1) the M31 RGB DD051 probability distribution func- 
tion (PDFrgb, Section [221), convolved with the DDOSl 
selection function, and (2) the MW dwarf star DD051 
PDF (PDFdwarf) convolved with the DD051 selection 
function. The DDOSl PDFs for M31 and MW stars are 
shown in Figure [2][b) . 



DSF(/i 



DD051 



target (/dD051 J 



Cddo 



-^total(/DD05l) 

/PDFrgb DSF d/DDOsi 

/ PDFdwarf DSF C?/dD051 



(2) 
(3) 



The scaled M31 RGB counts in each field are divided by 
the DDOSl correction factor to account for differences in 
the DDOSl selection function between fields (Equation 
[1]). For fields in which DDOSl photometry was not used 
for designing the masks, Cddo is set to unity. 

In field aO (i?proj ^ 30 kpc), an inner field with many 
high priority targets, the DDOSl-based selection of spec- 
troscopic targets increased the measured RGB / dwarf ra- 
tio by a factor of 3.4. In the outermost fields where 
there are only a few bright stars (/q < 23.0) with high 
DDOSl parameters per mask, the DDOSl-based selec- 
tion increased the measured RGB/dwarf star ratio by a 
factor of 1.4. 

The application of the DDOSl correction factor to 
the M31 RGB/MW dwarf ratio results in excellent 
agreement between the surface brightness estimates of 
fields with and without DDOSl-based target selection 



at R-o 



30 kpc (e.g., 'mask4' (no DDOSl) and 'aO' 



and 'a3' (DDOSl)). Interior to this radius fields did not 
have DDOSl photometry, while exterior to this radius 
the majority of the fields do have DDOSl photometry. 



The DDOSl correction factor removes the need to empir- 
ically adjust the normalizations of the surface brightness 
estimates for fields with and without DDOSl photome- 
try. 

3.2.2. Sampling of the Stellar Luminosity Function 

In designing the spectroscopic masks, higher priority 
is given to brighter targets because recovering a velocity 
for these objects is more likely. However, a mask at large 
-Rproj typically contains a higher fraction of faint targets 
than a mask at small i?proj due to the sparseness of stars 
in M31's outer halo. For example, the number of 20. S < 
/o < 22 targets is about twice that of 22 < /q < 22. S 
targets in the i?proj = 30 and 60 kpc fields aO and al3, 
while these two magnitude ranges contain comparable 
numbers of targets in the outer fields m8 and mil (at 
120 and 16S kpc). 

The M31 RGB /-band LF rises more steeply toward 
faint magnitudes than the MW dwarf /-band LF. There- 
fore the measured RGB /dwarf star ratio is expected to be 
higher in fields with a larger number of faint targets. If 
uncorrected, the larger number of faint targets in the out- 
ermost fields would bias the measured surface brightness 
profile to be flatter than the true profile. The amount of 
bias is limited because the spectroscopic targets span a 
relatively small apparent magnitude range and the radial 
velocity measurement fails for many stars with /q > 22, 
due to low S/N in the continuum. 

In addition, variations in observing conditions (e.g., 
seeing and transparency) from one mask to another can 
lead to variations in the completeness function of radial 
velocity measurements. We attempted to minimize this 
effect by increasing the total exposure time for masks ob- 
served under sub-optimal conditions. Nevertheless there 
are field-to-field variations in the faint-end limiting mag- 
nitude of the spectra, which in turn leads to slight vari- 
ations in the fraction of the LF that is included in our 
statistics. 

To account for the above effects, we estimate the fac- 
tor by which the measured M31 RGB/MW dwarf star 
ratio is increased in each field due to (1) a larger fraction 
of faint targets and/or (2) a fainter limiting magnitude 
for recovery of radial velocities from the spectra. The 
/ magnitude recovery function (IRF) is defined as the 
ratio of the number of stars with successful radial veloc- 
ity measurements (A^vei) to the total number of available 
targets in the photometric catalog (iVtotai), as a function 
of / magnitude. The correction factor (clf) for each field 
is defined as the ratio of two integrals: (1) the M31 RGB 
LF convolved with the / magnitude recovery function, 
and (2) the MW dwarf star LF convolved with the / 
magnitude recovery function. 



IRF(/) ^ 



_ J LFm31 IRF dl 
""^^ ~ J LFmw IRF dl 



(4) 
(5) 



This is analogous to the procedure used for calculat- 
ing the DDOSl selection efficiency factor (Section l3.2.1|) . 
The LF of M31 RGB stars is assumed to be dN/dl = 
j^qO.s/o^ with a cutoff at Iq = 20. S corresponding to the 
tip of the RGB. The MW dwarf star LF is assumed to 
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be a constantly rising lunction over the magnitude range 
of the spectroscopic survey; the slo pe is based on the B e- 
sangon Galactic population model ([Robin et al.ll2003l ). 

The scaled M31 RGB counts in each field are divided 
by the correction factor clf to account for differences in 
the target selection and limiting magnitude of the spec- 
troscopic observations in each field. The correction for 
sampling of the LF is a smaller correction than the cor- 
rection for DD051 selection: while the maximum vari- 
ation in the DD051 correction is a factor of ~ 4, the 
maximum variation in the LF correction is a factor of 
- 2. 

3.3. Other Sources of Error 

There are additional potential sources of error in the 
surface brightness estimates. We briefly discuss these 
and their possible impacts on the results below. In gen- 
eral, their effect on the measured surface brightness pro- 
file will be to add an extra source of noise in the mea- 
surements. 

3.3.1. Masks Targeting dSphs 

Some of the M31 fields target dSph galaxies in M31's 
halo (Section l2.3|) . In the dSph masks, the precision with 
which the number of field stars in M31's halo can be 
estimated depends on both the spatial extent of the dSph 
and the tightness of the locus of dSph member stars in 
the wios " [Fe/H] plane. The amount of parameter space 
occupied by dSph members limits the parameter space 
where M31 field stars can be securely identified. In most 
cases, the parameter space spanned by the dSph is small 
compared to the parameter space spanned by M31 halo 
stars. The magnitude of the error in M31 RGB star 
counts thus introduced can be estimated by counting the 
number of M31 stars in the range of [Fe/H] values and 
width in velocity occupied by dSph members, but on 
the opposite side of M31's Gaussian velocity distribution 
(i.e., if the dSph has a mean velocity of —450 km s^^, 
one can count the number of M31 stars at —150 km s~^ 
that would have been within the area of the v\os - [Fe/H] 
locus used to identify dSph stars). This additional source 
of error is generally found to be significantly smaller than 
the Poisson error; at the most, it is comparable in fields 
with very few M31 RGB stars. 

3.3.2. Population Gradients 

The PDFs used to identify M31 RGB stars were de- 
fined using a training set that was primarily drawn from 
fields in the inner, metal-rich regions of M31's stellar 
halo. If the RGB classification method is significantly bi- 
ased against selecting metal-poor stars, it is possible that 
large scale radial gradients in the metallicity of the halo 
population might bias the profile measurement. M31's 
stellar halo does grow increas i ngly metal-poor with in- 
creasmg radi us (iKahrai et all [2006at iKoch et al.l 120081: 
iTanaka et al.l 120101 1 A gradient in the age of the stel- 
lar population would also produce a bias, as changes in 
age produce similar color offsets as changes in metallicity 
for stars on the RGB. 

However, there does not appear to be a significant bias 
against identifying metal-poor RGB stars as M31 stars 
in the sample. If the color-based and metallicity-based 
diagnostics biased the selection against blue, metal-poor 



stars, it would be most apparent in the velocity range 
-200 to -150 km 8"^ This is where the PDFs of the 
M31 RGB and MW dwarf velocity diagnostics overlap 
(Figure [2]) and the effect of the velocity diagnostic on the 
overall likelihood is minimized. Figure [3] shows that the 
range of velocities spanned by stars classified as secure 
M31 red giants is similar whether the stars have [Fe/H]> 
-1 or [Fe/H]< -1. 

3.3.3. Photometric Data Quality 

Field-to-field variations in photometric accuracy and 
image quality lead to slight variations in the selection 
efficiency of M31 RGB spectroscopic candidates and the 
fraction of contaminating foreground MW dwarf stars 
and compact background galaxies observed on the spec- 
troscopic masks. In general, the limiting magnitude of 
the photometric data is significantly fainter than that of 
the spectroscopic data, thus the error introduced by this 
effect is expected to be small compared to the Poisson 
errors. 

3.3.4. Galactic Star-count Model 

The surface density of MW dwarf stars predicted by 
the B esancon Galactic population model (jRobin et al.l 
l2003f ) only enters the final surface brightness estimates 
in a relative sense. It is used solely as a normalization 
factor (A'bgm) to account for the changing density of 
MW dwarf stars across the survey, and a separate nor- 
malization factor is used to convert the surface densities 
to surface brightnesses (Section |4|). Therefore, an abso- 
lute error in the total number of stars predicted by the 
Besangon model would not affect the resulting surface 
brightness estimates. 

However, the Besangon model provides only a first 
order approximation for the true changes in MW stel- 
lar density across the footprint of the survey. A field- 
dependent error in the predicted projected density of 
MW dwarf stars could affect the resulting surface bright- 
ness profile. 

If the Galactic latitude dependence of the number 
counts in the model is incorrect, this would result in a 
systematic error in the surface brightness estimates with 
a magnitude that varies as a function of the Galactic lat- 
itude of the field. A factor of two systematic error in the 
Besangon model between the highest and lowest Galac- 
tic latitude fields would result in a maximum system- 
atic error of 0.75 mag; a factor of 1.5 error in the model 
counts would result in a maximum systematic error of 
0.4 mag, which is comparable to the Poisson uncertain- 
ties for many of the fields. 

However, there is not a direct correspondence between 
the projected distance of the fields from M31's center 
and their Galactic latitude. Thus, the primary effect of 
a systematic error in the Galactic latitude dependence 
of the Besangon model number counts will be increased 
scatter in the surface brightness profile rather than a bias 
in the slope of the profile fit. 

The Besangon model accounts only for the smooth 
components of the MW. However, the MW's halo is 
known to have abundant substructure, and there is one 
known substructure, the T riangulum- Andromeda feature 
()Rocha-Pinto et all I2004D , which is in the direction of 
M31 and may affect some of the fields. Since the in- 
creased MW stellar density due to substructure is not 
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accounted for by A^bgm, the subsequent increase in ob- 
served MW star counts would result in a fainter M31 halo 
surface brightness estimate for a line of sight that passes 
through MW substructure. Therefore, the net effect of 
significant MW halo substructure will be increased scat- 
ter in the observed M31 halo surface brightness profile. 

4. SURFACE BRIGHTNESS PROFILE OF MBl'S 
STELLAR HALO 

Figure [B] displays the surface brightness estimates (Ta- 
ble [21) as a function of projected distance from M31's 
center for each of the M31 fields shown in Figure [T] The 
surface brightness profile of M31's stellar halo clearly ex- 
tends to large projected distances (i?proj <; 175 kpc), with 
no evidence of a downward break in the profile (i.e., a 
transition to a steeper profile) at large radii. 

To take into account the flattened proflle of MSl's in- 
ner regions, the i?proj values of a small subset of the fields 
were converted to an effective radial distance along M31's 
minor axis. This effective i?proj was calculated adopting 
a flattening of h/a = 0.6 for i ?r>rni < 30 kpc, as mea- 
surcd for M31's inner spheroid ("Wal terbos fc Kennicuttl 
[1988 : Pritchet & van den Bcrgh 199^7" No correction 
was necessary for the majority of the innermost flelds 
(-Rproj ^ 30 kpc), since most lie along M31's minor axis. 
There are no strong observational constraints on the el- 
lipticity of MSl's spheroid at large distances, therefore 
the spheroid of M31 is assumed to be spherical {h/a — 1) 
beyond _Rproj = 30 kpc. This assumption is tested in 
Section [131 

Approximately a third of the fields targeted dur- 
ing the course of the SPLASH Keck/DEIMOS survey 
have c onclusive kinematical evidence of tidal debris (Ta- 
ble [Tl iGu hathakurta et all 120061 : iKalirai et all l2006bl : 
IGilbert et al. 2007, 200^^ In many cases, these kine- 
matically cold components can be associated with known 
tidal debris features identified as photometric overdensi- 
ties in star count maps. The left panel of Figure HI shows 
the surface brightness estimates including all M31 halo 
stars, while the right panel shows the surface bright- 
ness estimates with the kinematically cold tidal debris 
features statistically subtracted (Section l3.ip . In fields 
with kinematically cold tidal debris features, the error 
bars shown in Figure [S^b) include the uncertainty in the 
number of RGB stars belonging to the kinematically hot 
component, estimated from the maximum-likelihood fits 
to the field's stellar velocity distribution (Section [31). 

Statistical subtraction of tidal debris features reduces 
the scatter about the best-fit profile for the inner and 
medial fields; this effect is quantified below (Section l4.ip . 
Significant scatter in the surface brightness estimates is 
noticeable at i?proj ^ 90 kpc (Figure [SI), particularly in 
the right panel where tidal debris features have been sta- 
tistically subtracted from the M31 RGB counts. The 
most distant field in which kinematically cold compo- 
nents can be identified in the spectroscopic data is at 
-Rproj = 85 kpc (Table [H) . In fields at these large pro- 
jected distances, there are typically < 10 secure M31 
RGB stars per field. This makes discovering a tidal de- 
bris feature in the spectroscopic data prohibitively diffi- 
cult unless it happens to dominate the line of sight. It 
is likely that there is substructure in some of these outer 
fields that cannot be identified as such due to the lim- 
ited number of M31 stars. For example, fields 'streamE' 



and 'streamF' were placed on phot ometric overdensitie s 
identified as tidal debris features bv lTanaka et al.l(|2010D . 
However there is no definitive evidence of kinematically 
cold components in the spectroscopic data. 

In the following subsections, we will discuss fits to 
the surface brightness profile of M31 with and without 
the inclusion of kinematically-identified tidal debris (Sec- 
tion 14. ip and explore the extent of MSl's stellar halo 
(Section 14.21) including an examination of the M31 RGB 
sample for signs of MW dwarf star contamination (Sec- 
tion H^IIl). We also discuss the scatter of the data about 
the profile fit (Section 14. Sp . measure the ellipticity of 
MSl's stellar halo (Section 14. 4p . and compare our sur- 
face brightness estimates with M31 halo profile fits from 
the literature (Section l4?5l) . 

4.1. Fits to Profiles With and Without Tidal Debris 
Features 

The innermost data points (i?proj^ 20 kpc) are well 
within the region of MSl's halo where previous studies 
have found the profile to be consistent with an exten- 
sion of MSl's bulge (Section [1]). Therefore, we do not 
assume that contributions to the surface brightness from 
MSl's bulge and disk are non-negligible in these fields. A 
combined bulge, disk and power-law halo profile is fit to 
the MSI surface brightness estimates (Section [31) using a 
maximum-likelihood analysis. The intensity at a given 
radius is given by the functional form 

/ ^ /b(i?p..oj) + /d(i^proj) +Io( . !f ^2 ) ° (6) 

VI + [Kproi/ahr J 

where h and Id represent the profile of MSl's bulge 
and disk along the minor axis, —2a is the halo's two- 
dimensional power-law index, ah is a core radius inside 
of which the halo profile flattens, and Iq and are nor- 
malization constants. 

The recently pub l ished MSI surface brightness profile 
by iCourteau et al.l (|201l F) is used to fix the profi les of 
MSl's bulge (/b) and disk (/d). ICourteau et al.l ana- 
lyzed photometric observations of MSl's inner regions, 
and compiled them with previously published surface 
brightness estimates to measure MSl's surface bright- 
ness profile along the major and minor axes. Here, we 
use their preferred minor-axis fit "U," which included 
a Sersic bulge (n — 1.9, Be = 0.73 kpc), exponential 
disk (with scale length = 5.0 kpc), and power-law 
halo. While this particular bulge/disk decomposition re- 
produces the observed surface brightness profile along 
the minor axis, the kinematical data in the innermost 
fields do not sup port the pr esence of such a prominent 
disk component ([Gilbert et a l. 2007); an attempt to in- 
clude the kinematical data as a constraint on the inner 
minor-axis profile will be the subject of a future paper 
(Dorman et al., in preparation). 

Althoug h this work has s ignificantly more data at large 
radii than Courteau et al.l . their minor axis profile has a 
much higher density of points out to 30 kpc. Their pro- 
file thus provides significantly more leverage for fitting 
the transition from MSl's inner regions (bulge and disk) 
to MSl's halo than the data presented here. Therefore, 
we adopt the co re radius of a;i = 5.2 kpc found for the 
ICourteau et al.l 'U' fit. This parameter sets the ampli- 
tude of the halo at small radii, and affects the radius at 
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Fig. 6. — Surface brightness profile of M31's halo (Section |4]|, based on counts of spectroscopically confirmed RGB stars (Section [Sj . 
(a) Surface brightness estimates based on all M31 RGB stars in a given field, (b) Surface brightness estimates of the kinematically 
hot population of M31 RGB stars in each field (Section [sTTJl; due to the sparse nature of the outermost regions, tidal debris can only 
be identified in fields within Hproj $ 90 kpc. The statistical subtraction of identified tidal debris features results in a smoother surface 
brightness profile for MSl's stellar halo. Error bars are the combined Poisson errors of the M31 RGB and MW dwarf star counts. For 
fields with substructure statistically removed, the error in the M31 RGB counts includes the uncertainty in the fraction of the population 
belonging to the kinematically hot component. The solid curve shows the maximum-likelihood minor axis surface brightness profile fit to 
the data and includes a Sersic bulge, exponential disk, and power-law halo. The bul ge and disk com ponent pa ram eters are set to those 
found to best fit a compilation of observations of M31's light distribution by Courteau et al.l Il2011j ) (Section 14.111 . Red points denote 
measurements from fields outside the southern quadrant of MSl's halo, confirming that M31's stellar halo ext ends t o large projected 
distances (-Rproj > 150 kpc) and follows the same power-law decline in surface brightness in all directions (Section 14. 2.21 1. 



which the halo profile turns over. We adopt i?, — 30 kpc, 
which is safely in the regime where the halo dominates 
the surface brightness profile; this choice is arbitrary and 
only affects the value of lo {Ih{R*) — h)- 

The total model intensity is compared to the surface 
brightness estimates to determine the power-law parame- 
ters (normalization, Iq, and power- law index, —2a) that 
provide the maximum-likelihood fit to the data. The 
errors in the surface brightness estimates are used as 
weights when performing the maximum-likelihood anal- 
ysis. When calculating the fit to the data set with 
kinematically identified tidal debris statistically removed 
(Figure ID^b)), the fit was restricted to the radial region 
over which each field has sufficient numbers of M31 stars 
to identify kinematical substructure (i?proj < 90 kpc). 

Given the significant scatter in the data, especially 
at large radius, we explore the effect of individual sur- 
face brightness estimates and their errors on the best- 
fit power-law index using Monte Carlo tests. Each trial 
was randomly drawn from the observed sample, allow- 
ing redrawing of individual data points. For each data 
point drawn, the number of M31 RGB and MW dwarf 
stars used to determine the surface brightness estimate 
are adjusted to a random deviate drawn from a Pois- 
son distribution. The surface brightness estimates are 
then recalculated, and the maximum-likelihood power- 
law profile determined; 1000 trials were performed. 

The data are best fit by a power-law profile with in- 
dex — 2.2 ± 0.2 and normalization /io = — 2.51og/o = 



28.15 ± 0.15 mag arcsec"^. The quoted error estimate 
is the weighted standard deviation of the distribution of 
best-fit power-law indices from the Monte Carlo tests. 
This is fo rmally consisten t with the best-fit power-law 
found by iCourteau et al.l (index —2.5 ± 0.2 and fiQ = 
28.07±0.05 mag arcsec"^). The best fit power-law profile 
index for the sample with kinematically cold tidal debris 
features removed is also —2.2 ±0.3, fully consistent with 
the index found using all stars, but with a fainter halo 
normalization: /xq = 28.58 ± 0.14 mag arcsec"^. 

The measured halo profiles are not very sensitive to 
the adopted bulge and disk profiles. Even if the surface 
brightness is assumed to be dominated by the halo for 
the full range of the data, and only a power-law compo- 
nent is fit, the best-fit index is —2.4 ± 0.2, consistent 
with the power-law fit including disk and bulge com- 
ponents. The power-law index is also fairly insensitive 
to the adopted core radius. If a smaller core radius, 
Qh = 2.7 kpc, is adopted, consistent with an M31 halo 
measurement base d on blue horizontal branch stars as 
tracers of the halo ([Williams et al.l[2012t ) , the best- fit in- 
dex is -2.3 ± 0.2. 

The best-fit index is also not strongly dependent on 
the radial range used for fitting the data. If the ra- 
dial range is restricted to i?proj > 20 kpc, the best-fit 
power-law index is —2.1 ± 0.3 for the sample including 
all M31 stars; the best-fit power-law index is —1.9 ± 0.4 
for i?proj > 35 kpc. Although these best-fit values are 
shallower than the best-fit profile to all data points. 
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the power-law indices are all consistent within the er- 
rors. The best-fit index for the sample with kinemati- 
cal substructure statistically removed in the radial range 
20 < -Rproj < 90 kpc is also consistent: —2.0 ± 0.5. 

4.2. Extent of M31 's Halo 

4.2.1. Constraints on Contamination from MW Dwarf 

Stars 

M31 RGB stars are detected as far as 175 kpc from 
the center of M31. The most distant field included in 
this paper, d7 (i?proj= 230 kpc), has objects identified 
as marginal M31 stars, but no secure M31 RGB star de- 
tections (Section [2.2p . Given the large projected radial 
distances of the outermost fields, it is prudent to explore 
the effects that significant MW dwarf star contamination 
would have on the surface brightness estimates. Fore- 
ground MW dwarf star contamination consists of two 
physical groups of stars: those belonging to the MWs 
disk and halo. 

Dwarf stars in the disk of the MW have heliocentric 
line-of-sight velocities near ^ km s~^ (Figure [S]). The 
surface density of MW disk stars increases as the galactic 
latitude of the field approaches 0°; this effect is signifi- 
cant over the area spanned by our spectroscopic survey 
(Figure [T|). If MW disk stars are a significant source of 
contamination in the M31 RGB sample, fields closer to 
the MW disk (at lower absolute galactic latitude) should 
have systematically higher surface brightness estimates 
than fields farther removed from the MW disk. No clear 
trend is seen in the data (Figure [71 left panel) , although 
each radial bin encompasses a rather large range of radii 
and hence intrinsic surface brightnesses. 

The effect of the variation of surface brightness as a 
function of radius can be accounted for by subtracting 
the measured /i/ from the best fit profile (Figure [71 right 
panel). Intrinsic scatter of the M31 halo's surface bright- 
ness about the profile fit will be independent of galactic 
latitude, producing a pure scatter plot. If significant 
contamination from MW disk stars was present in the 
sample, there would be a clear trend in fiat — Mobs as a 
function of Galactic latitude, with preferentially larger 
values at lower absolute Ga lactic latitudes. Ba sed on the 
Besangon Galactic model ([Robin et al.ll2003() . the den- 
sity of MW stars along the line of sight is expected to 
increase by a factor of 8 between the lowest and high- 
est galactic latitude fields. If an approximately constant 
fraction of the MW disk stars present in each field creeps 
into the M31 RGB sample, the expected difference in 
surface brightness would be 2.2 magnitudes. Restricting 
the range of galactic latitudes to where the majority of 
the data lie (—30 < & < —15) leads to an expected dif- 
ference of 1.3 magnitudes. There is no apparent trend 
in the data even in the two outermost radial bins, let 
alone a trend this large. Therefore, we conclude that 
MW disk star contamination is not significantly affect- 
ing the M31 RGB surface brightness estimates. This also 
further demonstrates that the method of separating M31 
RGB and MW dwarf stars (Section [221) is highly effec- 
tive. 

Bright, distant main sequence turn-off stars in the 
MWs halo have a large range of velocities, and they 
can have very negative (M31-like) line-of-sight velocities 
(Figure[3l). Their relatively blue colors make them partic- 



ularly difficult to distinguish from M31 RGB stars, since 
the diagnostics described in Section [Z!2l have less power at 
bluer colors (Figure [SJ . Foreground MW halo stars are 
expected to have approximately constant surface den- 
sity across the area spanned by our M31 spectroscopic 
survey. Significant contamination from foreground MW 
halo stars in the M31 RGB halo sample would therefore 
manifest itself as a relatively constant (additive) surface 
brightness term, irrespective of radius. Naively, a con- 
stant term would simply be adjusted for in the empirical 
normalization of the data. However, given the spectro- 
scopic selection function, we observe more MW halo stars 
on masks at large radii (because more filler targets are 
placed on the masks), and contamination is therefore ex- 
pected to be larger in the outermost regions. If the term 
were large, it would therefore decrease the slope of the 
observed surface brightness profile. 

At some large radius, this term would become larger 
than the signal of actual M31 RGB stars, and the 
measurements would become foreground-limited, which 
would manifest as a constant observed surface bright- 
ness with increasing radius. However, the surface bright- 
ness estimates continue to decrease with radius, consis- 
tent with the power-law profile (although with increas- 
ing scatter) , out to the largest radii probed in our survey 
(Figure ini). Furthermore, the M31 halo sample and the 
likely MW halo stars (stars significantly bluer than the 
most metal-poor RGB isochrone. Figure [Sj occupy dis- 
tinct regions of the [Fe/H] -wjos plane. 

Finally, a comparison of the kinematics of the M31 halo 
sample and the likely MW halo stars shows that the two 
have very different distributions at all radii: the M31 
halo sample remains centered at the systemic velocity of 
M31, vios = -300 km s~\ while the blue MW stars dis- 
play a broad range of velocities approximately centered 
at the local standard of rest, visr = —175 km s~^. This 
is shown quantitatively in Figure [51 which displays the 
surface brightness estimates of a subset of the M31 spec- 
troscopic fields, calculated after splitting the M31 RGB 
sample at the systemic velocity of M31 (—300 km s^^). 
Incompleteness in the sample, which is expected to occur 
at velocities ov erlapping the MW s tar velocity distribu- 
tion (Figure [51 [Gilbert et all I2007D will cause points to 
fall above the one-to-one line, as surface brightness esti- 
mates calculated using the number of stars with v > vm31 
will be fainter than estimates calculated using the num- 
ber of stars with v < wm31- Conversely, contamination 
from MW halo stars will cause the points to fall below 
the one-to-one line. Since the same velocity PDFs are 
applied regardless of radius, incompleteness is expected 
to be roughly uniform, and is implicitly accounted for in 
the empirical normalization. Contamination from MW 
halo stars, however, is likely to be larger at faint surface 
brightnesses, where the ratio of M31 RGB stars to MW 
stars is low and the M31 RGB stars have relatively blue 
colors (Table [H Figures [H and [3]) . Reassuringly, the data 
points in Figure [8l are consistent with the one-to-one line 
even at faint surface brightnesses. 

Combined, these tests indicate that the M31 RGB 
counts are not significantly contaminated by misidenti- 
fied MW stars, and that we have not yet reached the 
radius at which the surface brightness estimates become 
foreground- limited. 
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Fig. 7. — Surface brightness of the M31 spectroscopic fields as a function of Galactic latitude. The left panel (a) shows measured surface 
brightness estimates, while the right panel (b) shows the scatter of the surface brightness estimates about the profile fit. The symbols used 
vary according to the radial distance of the field from M31's center. The dotted line denotes the Galactic Latitude of M31's center. If 
contamination from foreground MW disk stars was a significant factor in the M31 RGB surface brightness measurements, fields at lower 
absolute galactic latitude (closer to the disk of the MW) would be expected to have systematically higher surface brightness estimates 
than fields at comparable radii but higher absolute galactic latitude. Significant contamination would be most noticeable in the fields at 
large projected distances in M31's halo, where the number of detected M31 RGB stars is small. No trend is seen in the data, indicating 
that for egroun d MW disk stars are not a significant source of contamination in the surface brightness estimates of the outer halo fields 
(Section lOH l. 



4.2.2. Beyond the Southern Quadrant 

Most previous studies of the surface-brightness pro- 
file of the outermost regions of M31's stellar halo have 
been confined to the southern quadrant or M31's south- 
east minor axis (iGuhathaku rta ct al.' '20 051 llrwin et all 
120051 llbata et all I2007D . GuhathakurtaeLaU UOQ^ 
presented a surface brightness profile extending out to 
-Rproj ~ 165 kpc, generated from observations in fields 
in M31's so uthern quadrant, p rimarily along MBl's SE 
minor axis, llrwin et al.l ()2005[ ) likewise presented data 
along MBl's SE minor axis to pr ojected distances of 
-Rproj — 55 kpc. llbata et all (|2007D presented a photo- 
metric survey of MSl's entire southern quadrant out to 
-Rproj 150 kpc and detected substructure throughout a 
large portion of this region, much of it related to MSl's 
giant southern stream. These surface brightness profiles 
were all consistent with an power-law stellar halo. 

The stars in field mil, located on MSl's SE minor axis 
at a projected radial distance of ~ 165 kpc from MSl's 
center, are at a projected radial distance of ~ 50 kpc 
from the center of MSS. Although it is more plausible 
that the stars in field mil gravitationally belong to MSI 
given the relative sizes of the two galaxies, their velocity 
distribution do es not rule out the p ossibility that they 
belong to MSS (jGuhathakurta et al.l [2005) . In addition, 
a large H I bridge between M31 and M33 has been ob- 
serve d roughly along MSl's minor axis (|Braun fc Thilkeri 
120041) . and a possible interaction between MSI and MSS 
has been posited based on relatively high star count den- 
sities observed extending from MSS in the direction of 



MSI (iMcConnachie et al.l[2009l ). 

It could therefore be argued that some or all of the out- 
ermost halo MSI RGB stars (i?pioj ^ 100 kpc) in fields in 
the southern quadrant may be part of large extended sub- 
structures connected either with the MSl/MSS H I bridge 
or the giant southern stream, and not part of a global 
MSI halo extending to at least 170 kpc. If this were the 
case, significant enhancements would be expected in the 
surface brightnesses of fields on the SE minor axis and 
in the S quadrant compared to fields in one of the other 
quadrants of MSl's stellar halo. The data do not show 
any evidence of such an overdens i ty (Fi gure |6l). 

PAndAS (IMcConnachie et al.l I2009D has used the 
CFHT and the MegaCam instrument to survey the 
eastern, northern, and western quadrants of MSl's 
halo out to a ni a ximum distance of ^ 150 kpc. 
IMcConnachie et"aLl ()2009() presented the first evidence 
that the southern quadrant of MSl's stellar halo, al- 
though containing large amounts of substructure, still 
typified the extent and global structure of the full halo. 
Although they did not measure a surface brightness pro- 
file, they presented a stellar density map showing stars 
consistent with being MSI RGB stars from the SE mi- 
nor axis counter-clockwise to the NW minor axis. They 
also found significant photometric substructure over the 
full survey area. Tanaka et al. (2010) presented the first 
surface brightness profile composed of stars beyond the 
southern quadrant, using Subaru/Suprime-Cam photom- 
etry of contiguous fields along the NW and SE minor axes 



out to R 



proj 



100 kpc, and found reasonable consistency 



between the two profiles (Table [3]). 
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Fig. 8. — Surface brightness estimates of a subset of the M31 
spectroscopic fields, calculated after splitting the M31 RGB sam- 
ple at the systemic velocity of M31 (i>m31 = —300 km s"-'). 
Only fields with a reasonable expectation of symmetry in the 
M31 velocity distribution are shown here; this excludes fields with 
kinematically identified tidal debris features and containing dwarf 
spheroidal galaxies (Table [!}• The dotted line denotes the one-to- 
one line. MW halo stars have a broad range of velocities centered 
at 11 ~ —175 km s~^. Thus, significant contamination from MW 
halo stars would result in surface brightness estimates calculated 
using the number of stars with v > i;m31 to be brighter than esti- 
mates calculated using the number of stars with v < i'm31i and the 
points would fall below the one-to-one line. However, the points 
are consistent with the one-to-one line, indicating that the sample 
is not significantly co ntamin ated by MW halo stars at faint surface 
brightnesses (Section 14. 2. 11 1. 

The profile presented in Figure [6] confirms that M31's 
stellar halo follows a consistent power-law profile regard- 
less of position angle. Eleven of the M31 halo fields 
are on the northern side of the galaxy, spanning a large 
range of position angle and distance from MSl's center 
(Figure [1]). The surface brightness estimates from these 
fields are shown in red in Figure [6] and are fully consis- 
tent with estimates from fields along the SE minor axis 
and in the southern quadrant. Furthermore, secure M31 
RGB stars are found in both northern i?proj ~ 170 kpc 
fields. This supports the interpretation that the dynam- 
ically hot component in the southern quadrant fields are 
part of a global halo population that extends to at least 
170 kpc from MBl's center. 

4.3. Scatter of the Data about the Profile Fit: Effect of 
Substructure 

Although the best-fit power-law indices are consistent 
for both surface-brightness profiles shown in Figure [6l 
the surface brightness estimates including all M31 halo 
stars show slightly more scatter at any given radius 
(^proj^ 90 kpc) than the estimates that include only 
the fraction of stars belonging to the dynamically hot 
component of M31's stellar halo. This is quantified in 
Figures m^a) and (b). Statistical removal of stars associ- 
ated with kinematically cold tidal debris features results 
in a decrease of a factor of 1.3 in the root- mean-square de- 
viation of the data points about the best-fit profile. The 
rms value is computed using only data at -Rpioj < 90 kpc 
since that is the radial limit at which kinematical sub- 
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Fig. 9. — Comparison of the scatter of the surface brightness 
estimates about the profiles shown in Figure |6l as a function of 
distance from the center of M31 (loft panels) and in histogram 
form (right panels), (a) Values of (defined as /Xflt — fiobs) iov 
M31's surface brightness profile including all M31 IIGB stars (Fig- 
ure [Gja)). The histograms show the scatter about the best-fit af- 
ter controlling for the observational errors (A/^/cr^). (b) Values 
of Afi for the dynamically hot M31 RGB component and associ- 
ated profile fit (Figure |6l ft)); dynamically cold tidal debris features 
have been statistically removed, (c) The aggregate effect of sub- 
structure on the profile: values of A^ calculated with the surface 
brightness estimates including all M31 RGB stars (Figure[6la)) and 
the profile fit to the dynamically hot M31 RGB population (Fig- 
ure |5Ib)). Shaded histograms are for fields within Bproj < 90 kpc, 
the outer limit at which we are able to identify substructure. In the 
left panels, the dotted lines show the mean and root-mean-square 
deviation from the mean of the Afi values (blue: calculated for 
Kproj < 90 kpc, red: calculated using all surface brightness esti- 
mates). Statistical removal of stars associated with dynamically 
cold tidal debris features in M31's stellar halo results in a factor 
of 1.3 decrease in the rms devi ation of the data from the best-fit 
profile (panels a and b, Section l4.3l l. 

structure can be detected in the current data set. 

The mean enhancement of the surface brightness in 
a field due to kinematically identified substructure is 
1 mag; the root-mean-square deviation of the surface 
brightness enhancement in fields with kinematically iden- 
tified substructure is 0.5 magnitudes. Figure [9ljc) shows 
the aggregate effect of the inclusion of tidal debris fea- 
tures on the profile; it compares the surface brightness 
estimates including all stars with the best-fit profile to 
the surface brightness estimates with tidal debris statis- 
tically subtracted. 

Next, we investigate whether the scatter of the sur- 
face brightness of individual sight lines in M31's stellar 
halo increases with distance from M31. Physically, an 
increase in scatter with increasing i?proj could be caused 
by an increasing fraction of stars in the outer halo be- 
longing to substructure. Observationally, an increase in 
scatter with -Rproj could be due to the smaller number of 
stars observed in the outer fields: this simultaneously in- 
creases the error in the surface brightness measurements 
and reduces our ability to identify and remove substruc- 
ture. 

We can control for increasing error in the surface 
brightness estimates with increasing i?proj by analyzing 
the distribution of Ayit/cr^ values. If the scatter about the 
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Fig. 10. — Surface brightness profile of M31's stellar halo; data 
points are color-coded by the field's position angle from M31's mi- 
nor axis. Kinematically identified tidal debris features have been 
statistically subtracted from the surface brightness estimates in 
fields at i?proj< 90 kpc. Open points use iJproj ior the abscissa, 
while the abscissa value of the smaller filled points (connected by 
arrows) is their projection onto the minor axis assuming a flat- 
tening of b/a = .6 in M31's inner spheroid (Rproi < 30 kpc; 
Secti on[4l Waltcrbos & Kennicutt 1988; Pritchct & van den Bcrgh 
119941) . There is no apparent systematic offset of the surface bright- 
ness estimates as a function of angular position, indicating that the 
data are cons isten t with circular isophotes for MSl's outer stellar 
halo (Section rot . 

best fit profile is due solely to the observational errors, 
this distribution should have a standard deviation of or- 
der unity. The standard deviation is computed in three 
large radial bins: i?proj < 30 kpc, 30 < i?proj < 90 kpc, 
and i?proj > 90 kpc. The standard deviation is ~ 2 for 
the fit to fields including substructure, regardless of ra- 
dius. For the surface brightness estimates with tidal de- 
bris features removed, the standard deviation is 1.5 in 
both the -Rproj < 30 kpc and 30 < i?proj < 90 kpc bins, 
and 1.7 for fields at i?proj > 90 kpc, where the number 
of stars per field is insufficient to identify kinematically 
cold components. 

There are several sources of observational error not 
captured in the formal error estimates that could increase 
the scatter in the surface brightness estimates. The frac- 
tion of stars in tidal debris features is determined by 
fitting multiple Gaussian components to a field's veloc- 
ity distributions (Section 13.11) . However, Gaussians are 
likely an imperfect approximation of the true velocity 
distribution of the tidal debris: if there are more stars 
in the tails of the velocity distribution of tidal debris 
features than in the modeled Gaussians, this could lead 
to residual scatter about the profile. In addition, sub- 
structure in the MW's halo or a large scale error in the 
accuracy of the Besangon model, which is used as a nor- 
malization factor, could also increase the scatter in the 
surface brightness estimates (Section 13. 3|) . 

4.4. Shape of M31 's Halo 

The surface brightness estimates cover a large range of 
projected distances and position angles in M31's stellar 
halo. If M31's stellar halo has flattened isophotes with 



the same orientation as M31's disk, data points at in- 
creasing position angles from the minor axis would have 
systematically higher surface brightness estimates than 
points at similar radii but near M31's minor axis. If 
M31 had flattened isophotes with an orientation perpen- 
dicular to that of M31's disk, the data would show sys- 
tematic offsets in the opposite sense: data points with 
large position angles from the disk minor axis would have 
systematically fainter surface brightness estimates than 
data points at similar radii and near M31's minor axis. 

We investigate the data for evidence of such an effect in 
Figure flOl The surface brightness estimates and best-flt 
proflle are for the sample with tidal debris features sta- 
tistically removed to prevent bias of the measurement by 
strong, large-scale substructure (such as the giant south- 
ern stream). The surface brightness estimates are color- 
coded by their position angle from the minor axis (de- 
flned such that 9 < 90°). There is no signiflcant, system- 
atic offset of any subset of the data points from the proflle 
fit, nor does there appear to be an increasing systematic 
offset with increasing position angle from the minor axis. 
The data appear to be consistent with a spherical halo 
in M31. 

We quantify the degree of ellipticity of M31's stellar 
halo via a maximum-likelihood analysis, under the con- 
straints that the shape of the halo is assumed to be con- 
stant with radius and that the major axis of M31's halo 
is aligned with either the major or the minor axis of 
M31's disk. Each data point was projected along M31's 
(disk) minor axis, assuming an axis ratio for the halo 
isophotes. The halo axis ratio is defined with respect to 
the M31 disk's minor and major axes: a halo axis ra- 
tio smaller than one indicates the halo's major axis is 
aligned with the disk major axis, while a halo axis ra- 
tio greater than one indicates the halo's major axis is 
aligned with the disk minor axis. The best-fit profile was 
calculated as described in Section 14.11 The maximum- 
likelihood halo axis ratio was determined by comparing 
the values of each fit. Since the analysis is performed 
on the sample with kinematically identified substructure 
statistically removed, only fields with initial values of 
^proj< 90 kpc are included in the profile fits and 
measurement; beyond this radius individual fields have 
too few stars to identify kinematical substructure. 

The maximum-likelihood axis ratio for M31's halo is 
LOG^Qg"'^. This corresponds to elliptical isophotes with 
b/a = 0.94, with the major axis of the halo aligned 
along the minor axis of M31's disk, consistent with a pro- 
late halo. Circular isophotes (consistent with a spherical 
halo) are within the 90% confidence limits, which en- 
compass axis ratios of 0.96 to 1.17. The 99% confidence 
limits span a range of halo axis ratios of 0.92 and 1.25. 
The best-fit power-law index for a halo axis ratio of 1.06 
is —2.2, the same as the power-law index of the nominal 
fit (Figured Section SH]). 

Finally, based on previous 

work (iWalterbos fc KennicuttI :1988t 

iPritchet fc van den Berg hi I1994D . we assumed an 
ellipticity in the surface brightness isophotes for fields 
with Rp^o}< 30 kpc for the nominal fits (Section |4]). In 
practice, this significantly affected only two of the data 
points. The reader can see the effect of this assumption 
in Figure 1101 the open points show the surface bright- 
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ness estimates using the field's i?proj, while the smaller 
filled points (connected with arrows to the -Rproj points) 
show the same surface brightness estimates, but using 
the effective radius along the minor axis assuming an 
isophotal flattening oi b/a ~ 0.6 (with the major axis 
aligned with the major axis of M31's disk). Based on the 
results of the above analysis, we also tested the profile 
fits using i?proj for all data points, rather than assuming 
flattened isophotes in the inner regions, and found fully 
consistent fits with those quoted in Section 14.11 It is 
possible that the inner regions of M31's stellar halo have 
a different ellipticity than the outer regions; however, 
given the limited number of data points at small radii, 
most of which are near the minor axis, we are unable to 
constrain the shape of MSl's inner halo separately from 
the shape of the outer halo. 

4.5. Comparison with Previously Published Profiles 

In the last several years multiple measurements of the 
surface brightness profile of MSl's outer halo have been 
made, primarily utilizing photometric observations along 
the southeast minor axis. These are summarized in Ta- 
ble 121 Although these studies largely use independent 
data sets and analysis, there is partial overlap in the 
footprints of several of the surveys; for example, seven 
of t he 39 spectros c opic fi elds presented here overla p with 
the iTanaka et all (|2010D footpri nt. The fCourtca u et al.l 
(|201l[ ) analvsis makes use of the 'Irwin et al. (2005^ data 
set and an updated Guhathakurta et al. (2005) data set. 

Figure [TT] shows the surface brightness estimates in log- 
linear form along with a range of literature fits to ob- 
servations of MSl's outer halo (-Rproj ^ 30 kpc) surface 
brightness profile, including power-law profiles, Hern- 
quist models (left panel), and exponential profiles (right 
panel). The Hernquist models provide reasonable fits to 
our data set; the smallest Hernqu ist model scale rad ius 
from the literature (Rs = 17.1 kpc. lTanaka et al.ll2010[ ) is 
the most consistent with our data set over its full radial 
range. The data are consistent with the exponential fits 
only within relatively small radial ranges. A power-law 
profile or Hernquist model clearly provides a better fit 
than an exponential profile over the full radial range of 
the data. 

The surface brightness estimates presented here ex- 
tend to larger projected distances and cover a larger 
azimuthal range in MSl's stellar halo than any of the 
previously published profiles. Our technique of using 
counts of spectroscopically confirmed stars to estimate 
a surface brightness has an advantage at large projected 
distances over estimates based on photometric observa- 
tions. In the sparsest outer regions of MSl's stellar halo, 
contaminants from the foreground (MW dwarf stars) and 
the background (distant unresolved galaxies) greatly out- 
number the target MSI halo population. Given the large 
area of sky covered by MSl's stellar halo (over 20° in 
diameter), and the strongly varying MW stellar density 
over this region of sky, the standard technique of iden- 
tifying a comparison field and statistically subtracting 
it from the science image to remove contamination docs 
not work. Spectroscopic observations identify individual 
objects as foreground MW stars, background galaxies, 
or M31 stars, enabling us to probe MSl's halo to large 
projected distances. 



5. DISCUSSION 

The density profiles found in si mulations of stellar 
halo formation via accre t ion fe.g.. ["Bullock fc JohnstonI 
[200l ICooper et alllMot iFont et al.l 12011 " are in good 
agreement with the meas ured surfac e brigh tness pro- 
file of MSl's stellar halo. iFont et al.l ()201lD presented 
a study of w 400 stellar halos, simulated using a cos- 
mological smoothed particle hydrodynamics code. Since 
stellar halo formation is a somewhat stochastic process, 
the exact properties of individual halos do differ. How- 
ever, the median spherically averag ed density profiles at 
large radii in the iFont et all ()201lD simulations follow a 
power-law profile with p(r) oc r~^'^ (leading to a surface 
density distribution of S(r) cx r~^'^), in good agreement 
with that m easured for the MSI halo s urface brightness 
profile. The iBullock fc JohnstonI ()2005D simulations pre- 
dicted stellar density profiles that closely followed Hern- 
quist models with scale radii in the range ~ 10 — 20 kpc, 
also in good agreement with the observed surface bright- 
ness profile of MSI fFigure lTlT a)). The published profiles 
typically show the aggregate density of the full stellar 
populations of the simulated stellar halos, including tidal 
streams, and are thus most comparable to the observed 
profile including all MSI stars. 

The data show no evidence for a downward break in 
the surface brightness profile at large radii. Most ACDM- 
based predictions for stellar ha l o profiles (IBullock et al.l 
I2OOII: IBullock fc Johnston|[2005l:ICooper et al.ll2010[ ) sug- 
gest that radially-averaged profiles of accreted stellar 
halos around MSl-size galaxies should start to break 
rapidly around r ~ 100 kpc, taking on a steeper pro- 
file than that interior to r ^ 100 kpc. Of the 1 1 sim- 
ulated stellar halos in IBullock fc JohnstonI (|2005| ) , only 
one, halo 9, did not break to a steeper density profile 
at large radius; it followed a Hernquist profile with scale 
radius of 18 kpc, which is sim ilar to the Hernquist scale 
radius for MSI measured by 'Tana ka et al.l (|20Tci) and 
shown to be consistent with the data in Figure [TT] This 
simulated halo had an unusually large fraction of recent 
low-mass accretion/disruption events. The large spread 
in surface brightness in the fields at large radii may indi- 
cate that MSI has undergone a similar recent accretion 
history. 

Simulations of stellar halo formation via accretion in 
a cosmological context also show that present-day stel- 
lar halos should contain ma ny photometrically i denti- 
fiable tidal streams (e.g., 'Bullock fc JohnstonI 120051 : 
John ston et"alll2008t iCooper et al. 20l3). Thus, a study 
utilizing data from isolated sight lines through a stellar 
halo, like this one, should see significant scatter in the 
surface brightness profile due to the presence or absence 
of tidal streams along each line-of-sight. Furthermore, 
the simulations find that the contrast of stars in coher- 
ent tidal debris streams to stars in the diffuse stellar halo 
should become greater at large radii, as these streams are 
typically due to more recent accretion events. The inner 
regions of the simulated stellar halos tend to be built 
from early, relatively major mergers, and are more likely 
to have undergone significant phase-mixing resulting in 
a smoother stellar distribution. 

The observations presented here provide one of the first 
direct observational tests of this aspect of the simula- 
tions. The surface brightness estimates show significant 
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Fig. 11. — Comparison of the surface brightness estimates of M31's stellar halo (as in Figurc|6la)) with a range of profile fits from the 
literature, in log-linear form. In addition to the best-fit power-law profile to the data (black solid curve; same a s in Figure [6ta)), alternative 
power-law indices and Hernquist models (left) and exponential fits (right) have been overlaid (Table|3l Section l4.5ll . (a) Alternative power- 
law indices for the outer halo component are overlaid on the data: —2.9 (red dotted curve) and —1.75 (blue dashed cur ve). A lso shown 
are two Hernquist model fits from the literature, (b) Exponential fits with scale lengths of 14 kpc (blue curve: Jrwin et ani2005l) . 31.6 kpc 
(red curve; Ibata et al. 2007), and 45.1 kpc (magenta curve; Ibata et al. 2007). The exponential fits are shown as solid lines over the radial 
range in which they were measured, and as dotted lines outside this radial range. Hernquist models are consistent with the full range of 
the data, while exponential profiles are only consistent with the data in limited radial ranges. 



ficld-to-ficld variation in the stellar content of M31's halo 
at all radii, beyond that expected based on the estimated 
observational errors. Tidal debris features are identi- 
fied via stellar kinematics in 1/2 of the fields within 
^proj ^ 90 kpc. The tidal debris features generally com- 
prise over half of the population in a given field, with the 
full range of measured substructure fractions ex tending 
from 2 9% in field d9 to 75% in field H13s. BullocklFii] 
(|200lD estimated the expected field-to-field variation in 
star counts in stellar halos built via the accretion of satel- 
lites to be a factor of ^ 2 to 3 at large radii, in rough 
agreement with the observed field-to-field variation in our 
data set. 

While removing kinematically identified substructure 
reduces the root-mean-square deviation of the data 
points about the profile fit, the amount of scatter in the 
data about the profile fit remains greater than that ex- 
pected from the estimated observational errors, suggest- 
ing that the outer halo is composed of many small ac- 
creted objects that are unresolved into kinematically cold 
clumps or photometrically distinct streams at the limit 
of the currently available data. Interestingly, the data do 
not show evidence of an increasing fraction of coherent 
tidal streams with radius, as expected from simulations 
and discussed above. It is worth noting, however, that 
the inner regions of M31's halo contain a high degree of 
substructure from the recent, relatively massive accre- 
tion event which formed the giant southern stream; in 
this sense, M31 may not be typical. 

A detailed comparison of the fraction of substructure 
with radius, its effect on the observed surface bright- 



ness profile, and an analysis of the field-to-field variation 
with that predicted in simulations should yield interest- 
ing constraints on the merger history of M31 and the 
properties of its former satellite populations (e.g., the 
LF) , and may help to provide constraints on the assump- 
tions for satellite properties used in simulations. This 
kind of detailed comparison is currently only possible in 
M31. M31's proximity allows us to obtain detailed ob- 
servations of the entire stellar halo, including the crucial 
kinematical data that allows us to separate tidal debris 
features from the smoother underlying stellar halo in in- 
dividual lines-of-sight. 

Interestingly, recent simulations employing smoothed 
particle hydrodynamics (and thus explicitly simulating 
the evolution of the baryons, including star formation) 
imply that a significant fraction of the stellar population 
in the inner regions of s tellar halos may be formed in 
the host galaxy potential (iZolotov et al.ll2009l : iFont et alJ 
120111: iMcCarthv et al.ll2012D . This could have several ef- 
fects on the observed stellar distribution, including nat- 
urally providing an inner spheroid (i?proj ^5 10-30 kpc) 
with different structural and chemical properties than 
the o uter spheroid, as has b een observed in M31. In- 
deed, IMcCarthv et al.l ()2012D showed that in their sim- 
ulations, the in situ stellar component results in a sig- 
nificantly flattened, oblate stellar density distribution 
in the inner regions (semi-minor to semi-major axis ra- 
tio of 6/a ~ 0.3) that slowly becomes more spheri- 
cal with increasing radius, reaching h/a ~ 1 at radii 
of 100 kpc. The inner regions (i?proj ^ 30 kpc) of 
M31's spheroid have been shown to have flattened surface 
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brigh t ness contours {b/a 0.6: iWalterbos fc Kennicutti 
[1981 iPritchet fc van den Berghi I1994D . while the outer 
regions appear consistent with circular surface bright- 
ness contours {b/a ~ 1.0; this work). In contrast, the 
TV-body plus semi-analytic models, in which star parti- 
cles are 'painted' onto the dark matter distribution, tend 
to fin d a prolate density distribution for the inner stellar 
halo (jCooper et al.|[2010l ). 

Detailed comparisons of the observed distribution of 
stars in M31's stellar halo with simulations of halo for- 
mation will be crucial for constraining the fraction of 
stars in stellar halos formed in situ. Spectroscopic stud- 
ies such as this one will have a natural advantage for such 
comparisons because the spatially diffuse, kinematically 
hot component can be separated from and compared to 
the tidal debris in every field. 

6. CONCLUSIONS 

We have analyzed data in 38 lines of sight in M31's 
stellar halo, including all quadrants and ranging from 
9 to 230 kpc in projected distance from M31's center. 
Spectroscopic observations of stars were obtained in a 
total of 108 multi-object slit masks. The stellar spectra 
were used to identify samples of M31 RGB stars and 
MW dwarf star contaminants, and to identify M31 stars 
belonging to tidal debris features. Counts of M31 RGB 
stars and MW dwarf stars were used to determine the 
surface brightness in each line-of-sight, yielding a surface 
brightness profile of M31's stellar halo. 

We expand upon our main conclusions below. In brief, 
our conclusions are as follows. 

• We present a clear detection of M31's stellar halo to 
projected distances of 175 kpc from M31's center. 

• Multiple tests indicate there is minimal contami- 
nation from misclassified foreground MW stars in 
the sample of M31 RGB stars. 

• The surface brightness profile follows a power-law 
with an index of —2.2. 

• The statistical subtraction of tidal debris allows us 
to investigate its effect on the measured surface 
brightness profile. 

We present a clear detection of MSI 's stellar halo 
to projected distances of 1 75 kpc from M31 's center 
(Figure\B^. Our data set includes lines-of-sight spread 
throughout all quadrants of M31's stellar halo, covering 
a range of radii. The data show that M31's stellar halo 
extends to at leas t ~ 2/3 of M31 's estimated virial ra- 
dius 260 kpc; S eigar et al.ll2008[ ). In aggregate, a large 
spread in the velocities of M31 halo stars is observed even 
in fields at large projected distances from M31's center 
(i?proj> 100 kpc, Figure 

Multiple diagnostic tests indicate there is minimal con- 
tamination from misclassified foreground MW stars in 
the sample of M31 RGB stars. The spectroscopic data 
enable us to remove foreground and background con- 
taminants, allowing us to trace M31's stellar halo to 
fainter surface brightnesses, and larger radii, than studies 
based on pho t ometr y alone. Tests show that using the 
I Gilbert et all ()2006D likelihood method to identify M31 
RGB and MW dwarf stars does an excellent job of cleanly 



separating the two populations (Section 14.2. 1[ Figures [7] 
and [8]). If the M31 RGB sample was significantly con- 
taminated by MW disk stars, there would be a strong 
dependence on the surface brightness of fields as a func- 
tion of Galactic latitude. No such dependence is seen 
in the data. The density of MW halo stars is expected 
to be relatively uniform across the face of M31, so sig- 
nificant contamination would lead to a constant surface 
brightness measurement. Our sample instead displays a 
systematic fall off with projected distance from M31's 
center. Furthermore, surface brightness estimates in in- 
dividual fields calculated by splitting the M31 sample at 
the systemic velocity of M31 are consistent even at faint 
surface brightnesses, indicating that contamination from 
MW halo stars is not significant. Even at large radii 
(^proj> 100 kpc) there are stars with properties consis- 
tent with M31 RGB stars rather than MW halo or disk 
stars (Figures [2] and 111) . 

The surface brightness profile follows a power law with 
an index of —2.2. This is consistent with what previous 
studies that have mainly targeted the southern quadrant 
of M31 have found. This work establishes that an ^ r~^ 
power- law profile is a good description of M31's outer 
halo (i?proj ^ 20-30 kpc) regardless of quadrant or ra- 
dius. The data show no evidence of a downward break 
in the profile at large radii (i?proj ^ 100 kpc). Since the 
fields cover a large range in both radii and position angle 
around M31, they provide leverage on the shape of M31's 
halo. The best-fit elliptical isophotes have b/a = 0.94, 
with the major axis of the halo aligned along the mi- 
nor axis of M31's disk, consistent with a prolate halo. 
However, a value of 6/a = 1 is well within the 90% con- 
fidence limits of the fit, so the data are consistent with a 
spherical halo. 

The statistical subtraction of tidal debris allows us to 
investigate its effect on the measured surface brightness 
profile. The ability to remove stars associated with tidal 
debris from individual lines of sight is a unique advantage 
of our spectroscopic survey over photometric surveys of 
M31's stellar halo. There are enough stars per field to 
discover kinematically cold tidal debris features for fields 
with -Rpioj ^ 90 kpc. Tidal debris features are identi- 
fied using kinematics in '--^ 1/3 of the spectroscopic fields 
(~ 1/2 of fields with i?proj ^ 90 kpc). In fields with sub- 
structure, the mean enhancement of the surface bright- 
ness is 1 mag, with a root-mean-square deviation of the 
sample about the best-fit profile of 0.5 mag. The scatter 
about the mean surface brightness profile decreases once 
kinematically identified tidal debris features are statisti- 
cally subtracted. This implies that there is a compara- 
tively diffuse stellar component to M31's stellar halo out 
to at least -Rpioj 90 kpc, with relatively prominent tidal 
debris features superimposed. 
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Basic Results of Spectroscopic Observations. 
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The distance of the field from M31's center, obtained by calculating the mean distance from M31 of all stellar spectra in a 
field. 

^ Bandpasses of the photometric catalog used to design the spectroscopic masks. 

Number of unique stellar spectra with successful velocity measurements, excluding duplicate measurements and alignment 
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"* These columns report the nu mber of stars designated as M31 RGB halo stars or MW dwarf stars as determined by the 
diagnostic method described in jl2.2l In dwarf galaxy fields, RGB stars identified as likely dSph members are not included in 
the M31 RGB count (Section [2^ 
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TABLE 2 

Surface Brightness Estimates Based on a 
Spectroscopically Confirmed Sample. 
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For fields in which only one mask was observed or all mask 
centers overlapped, the R.A. and decl. given is the center of the 
spectroscopic mask(s). For the rest of the fields, the R.A. and decl. 
are the mean R.A. and decl. of all observed masks. 
^ The normalization of t he surface brightness estimates is based 
on lCourteau et al.l I I2011I ) (Section [Hi . 

Empty values indicate that no kinematic substructure was de- 
tected in the field. 
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TABLE 3 

Recent Profile Measurements of M31's Halo. 



D t LIU J 


'-mm 


R 

J l-max 


P/^m/Of— 1 tnsr 


J_VCLU(X Lr y LH_^ 


Notes 

i. > U LCD 






rknr'l 
I, Kpc; 




(rn lypej 










(Scale Length)^ 






Guhathakurta et al. (2005) 


30 


165 


-2.3 


Spectroscopic 


S quadrant 












(included 10 of the fields presented here) 


Irwin et al. ('2005")'' 


20 


55 


-2.3 


Photometric 


SE minor axis 






(14) 


(Exponential) 




Ibata et al. f2007')'^ 


30-35 


90-130'' 


— 1.9 ± 0.12 


Photometric 


SE minor axis 








(31.6 ± 1.0) 


(Exponential) 










(54.6 ± 1.3) 


(Hernquist) 






90 


130 


(45.1 ±6.0) 


(Exponential) 










(53.1 ±3.5) 


(Hernquist) 




Tanaka et al. (2010")='" 


20 


100 


-1.75 ±0.13 


Photometric 


SE minor axis 








(22.4 ± 2.3) 


(Exponential) 










(31.7 ±6.7) 


(Hernquist) 






20 


100 


-2.17 ±0.15 




NW minor axis 








(18.8 ± 1.8) 


(Exponential) 










(17.1 ±4.7) 


(Hernquist) 




Courteau et al. (2011")* 





165 


-2.5 ±0.2 


Photometric and 


SE minor axis preferred fit 








-2 6+"'^ 
^•0-0.2 


Spectroscopic 


(included 11 of the fields presented here) 


Williams et al. f2012) 


2 


35 


Photometric 


BHB stars used as tracers of the halo 


This work 


9 


176 


-2.2 ±0.2 


Spectroscopic 


All quadrants, all M31 stars 




35 


176 


-1.9 ±0.4 




All quadrants, all M31 stars 




9 


90 


-2.2 ±0.3 




Tidal debris statistically subtracted 


^ The scale radii (in kpc) of 


exponential fits and Hernquist models 


are denoted within parentheses. 



^ The llrwin et al.l II2005I ) study contained integrated light measurements and RGB star count measurements of M31's surface 
brightness from M31's center out to 55 kpc. Only the data from 20 kpc outward was used to determine the profile of the outer 
halo. 



The profile fit is to metal-poor ([Fe/H] < —0.7) stars only. 

The profile is fit to stars in two radial regions. The data between 35 and 90 kpc were excluded from this analysis due to 
numerous tidal d ebris fe a tures. 

° Similar to Ibat a et al.l II2007I) . regions deemed overdense (i.e., strongly contaminated by tidal debris features) were excluded 
from the profile fit. 

* The [Courteau et al. ( 2011) study com bined integ rated light m easurements (Choi et al. 2002), photometry-based RGB star 
counts (from Pritchct & van den BcrgQ ( 119941 ) and llrwin et al. l (2005)), and spcctroscopy-bascd RGB star counts (11 of the 
fields presented here). They simultaneously fit for the bulge, disk, and outer halo profiles. Five of their data points were at 
^proj > 55 kpc. 



